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Nucleosynthesis in 7— ray bursts outflows 
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Abstract. It is shown that fusion of neutrons and protons to 4 He nuclei occurs in 7— ray burst outflows in a 
process similar to big-bang nucleosynthesis in the early Universe. Only the surviving free neutrons can then 
decouple kinematically from the charged fluid so that the mufti- GeV neutrino signal predicted from inelastic 
nuclear n— p collisions is significantly reduced. It is also argued that a sizeable fraction of ultra-high energy cosmic 
rays accelerated in 7— ray bursts should be 4 He nuclei. 
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1. Introduction. 

Multi-wavelength observations of 7— ray bursts in the past 
decade have given increasing evidence in favor of the "ex- 
panding fireball" model (Paczyhski 1986, 1990; Goodman 
1986, Shemi & Piran 1990; Piran 2000 and Meszaros 2002 
for reviews), in which a photon - pair plasma loaded with 
a small admixture of baryons expands relativistically and 
converts the initial energy into baryon kinetic energy. In 
the internal/external shock scenario (Rees & Meszaros 
1992; Rees & Meszaros 1994; Piran 2000 and Meszaros 
2002) this kinetic energy is dissipated in shocks, giving 
rise to the 7— ray burst phenomenon. The identity of the 
inner engine, the source of energy and the mechanism of 
collimation remain however unknown. 

One should not expect that only protons are injected 
in the accelerating wind. As a matter of fact, all theo- 
retical proposals for the progenitor involve compact ob- 
jects, e.g. neutron stars/black hole mergers, or imploding 
massive stars, in which the baryon load must be neutron 
rich. This has triggered recent interest in the study of phe- 
nomenological consequences of neutron loading in 7— ray 
bursts (e.g. Derishev et al. 1999a, b; Bahcall & Meszaros 
2000; Fuller et al. 2000; Meszaros & Rees 2000). 

One possible consequence is nucleosynthesis of neu- 
trons and protons to heavier nuclei. A recent study ar- 
gued that nucleosynthesis should not take place as the 
dynamical timescale (taken to be id yn ~ 10~ 5 sec) is 
too short (Derishev et al. 1999). However this particular 
value rather constitutes a strict lower limit to the dynam- 
ical timescale since it implies a source size < 3 km. This 
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timescale is nonetheless bounded from above by the short- 
est variability timescale observed w 10 -2 sec (Piran 1999). 

The purpose of this Letter is to show that nucleosyn- 
thesis to 4 He can occur and is generically efficient, pro- 
vided the dynamical timescale > 10 _4 sec (Section ||). 
Most available neutrons and protons are then bound into 
4 He and only the surviving neutrons can decouple kine- 
matically from the charged fluid component (Section ||). 
Thus this significantly weakens the 5 — 15 GeV neutrino 
signal expected from n — p inelastic collisions (Derishev 
et al. 1999a; Bahcall & Meszaros 2000; Meszaros & Rees 
2000). Furthermore, if 7— ray bursts accelerate particles to 
ultra-high energies > 10 20 eV (Levinson & Eichler 1993; 
Vietri 1995; Waxman 1995; Waxman 2001 for a review), 
the ultra-high energy cosmic ray spectrum should com- 
prise a sizeable fraction of 4 He nuclei (Section |h. 



2. Nucleosynthesis. 

The fireball wind can be modeled as a pair plasma with 
luminosity L = L50 x 10 50 ergs/s injected into a solid an- 
gle f2 = 4ir9 2 , where 9 2 = (l — cos #1/2) represents the 
effect of beaming into two jets of half-opening angle #1/2. 
Observations suggest a roughly uniform i 50 ~ 1 among 
7— ray bursts and a varying 9 with mean value 9 ~ 0.1 
(Frail et al. 2001; Panaitescu & Kumar 2001; Piran et al. 
2001); for a wind of typical duration ~ 10 sec, this yields a 
total equivalent isotropic output energy ~ 10 53 ergs. The 
bulk Lorentz factor is written T, and its saturation value 
77 = L/Mc 2 is the ratio of total luminosity to baryon out- 
flow (Meszaros et al. 1992). Injection takes place at ra- 
dius r a — r .7 x 10 7 cm, with initial temperature To ~ 
OmLl^g^r^MeV/kB, using 11/2 for the number 



2 



Martin Lemoine: Nucleosynthesis in 7— ray bursts outflows 



of degrees of freedom (photons + pairs). In the wind frame 
the cjccta is indeed similar to the early Universe before 
big-bang nucleosynthesis (Shemi & Piran 1990). However 
the physical conditions are very different: in a 7— ray burst 
ejecta, the comoving timescale ~ 10 _3 sec and the pho- 
ton to baryon ratio n 7 /nb — 4.1 x 10 4 L 5 ^ 4 O^^r^ 77300, 
with 7/300 = 7//300, to be compared with rjdyn ~ 100 sec 
and 77, 7 /?T,b ~ 10 10 for the early Universe. A large dynam- 
ical timescale favors nucleosynthesis but a high entropy 
acts against. Big-bang nucleosynthesis is also hampered 
by a small neutron to proton ratio due to neutron decay 
and late freeze-out of the weak interactions that intercon- 
vert protons and neutrons. Here one expects equal n,p 
mass fractions X n ~ X p if baryons come from photodis- 
sociated nuclei, and neutron decay is insignificant on a 
millisecond timescale. If T > 6.5MeV (id^/lO^sec)" 1 / 5 , 
the rate of weak interactions becomes larger than the 
fireball expansion rate and the neutron to proton ratio 
achieves equilibrium independently of the initial composi- 
tion, X n /X p — > exp[(77T, p — m n )c 2 /k^To\ ~ 1. However this 

occurs for < 0.02r^£l^ 2 (i dyn /10~ 3 sec) 2 / 5 , i.e. for the 
most highly beamed or compact 7— ray bursts, or those 
with the longest dynamical timescales (see also Fuller et 
al. 2000 for a discussion of the impact of a neutrino flux 
on the initial chemical composition of the outflow) . 

Adiabaticity during the acceleration phase implies that 
the bulk Lorentz factor T and the temperature T evolve 
as T oc r _1 oc r , with r the radial coordinate in the 
lab frame. If the flow can be approximated as quasi- 
one dimensional, the relativistic Bernoulli equation fur- 
ther implies r oc Or (Blandford & Rees 1974). In par- 
ticular Meszaros et al. (1993) use T w 0r/r o , which is 
obviously valid only in the limit r 0~ l r o . The exact 
numerical prefactor in this relation depends on the details 
of the injection phase at radius r ~ r a . In what follows, 
the early behavior will be phenomenologically modeled as 
r = (l + 0r/r o ), with T oc T -1 , which approaches the twin 
exhaust model of Blandford & Rees (1974), in which the 
flow is collimated by outside pressure forces and acceler- 
ates to relativistic velocity through a de Laval nozzle. A 
more exact solution requires solving the complex problem 
of injection and collimation in the early phase of the flow. 
Steps in this direction have been accomplished recently by 
Levinson & Eichler (2000) who studied the hydrodynamic 
collimation of a 7— ray burst outflow by a wind emanat- 
ing from a torus. Their solution indeed reproduces a nozzle 
with nearly constant cross-section at small radii similarly 
to the Blandford & Rees (1974) model. 

The dependence of the bulk Lorentz factor on r is 
important as it gives the dynamical expansion timescale 
idyn = r/Tc. Interestingly in this context the scaling 
r oc Or suggests that highly beamed outflows have a 
longer dynamical timescale, hence should not show varia- 
tions on short timescales due to erasure of inhomogeneities 
on scales smaller than the sound horizon < cidyn/%/3. 
Nevertheless, in order to circumvent uncertainties related 
to the modeling of the early evolution of the bulk Lorentz 
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Fig. 1. Mass fractions vs. radius for the case r Q = 3 x 
10 6 cm, = 0.1 (fi/47r = 10~ 2 ), dynamical timescale 
^dyn = 10~ 3 sec, L — 10 50 ergs.s -1 and initial mass com- 
position X n = 0.54, X p = 0.46. 



factor, all results that follow will be shown as a function 
of dynamical timescale tdyn instead of T. This is justi- 
fied as tdyn and the entropy are the two parameters that 
control the efficiency of nucleosynthesis. Moreover it was 
checked that a spherically symmetric wind with the sim- 
ple expansion law L — r/r a gives the same result than a 
jet with respect to nucleosynthesis provided the dynamical 
timescale and entropy are the same. 

In the ultra-relativistic regime T 3> 1 the tem- 
perature decreases exponentially fast with comoving 
time oc log(r/r ). In the outflow thermal equilibrium 
is ensured between all species all along nucleosynthe- 
sis, and neutrons remain coupled to protons through 
nuclear scattering with velocity averaged cross-section 
(<7 n _ p 7;) ~ 30mb c (Derishev et al. 1999b). Decoupling 
occurs when l/n p (a n - p v) > idyn, at temperature T w 
2i50 1/1 ^o / i 6 ^f 7 %30 / o( i dyn/10- 3 sec)- 1 /3 kc v/fc B , i.e., well 
after nucleosynthesis has taken place (see below). 

In order to estimate the final abundances of nuclei 
synthesized, use was made of a big-bang nucleosynthe- 
sis numerical code whose time and entropy evolutions 
were modified accordingly. This code accounts successfully 
for synthesis of elements up to 7 Be when compared to 
other big-bang nucleosynthesis calculations. An example 
of the outcome of nucleosynthesis is shown in Fig. |l| for a 
7— ray burst with fiducial parameters r a ^ = 0.3, = 0.1, 
L 50 = 1, initial composition X n — 0.54, X p = 0.46 
(photodissociated 56 Fe nuclei) and dynamical timescale 
tdyn — 10 -3 sec. The final mass fractions of elements pro- 
duced are X 4 ~ 0.84 (for 4 He), X n ~ 0.12, X p ~ 0.04, 
— 0.005 and other elements are produced only in 
traces. Although the numerical code used cannot deal with 
elements beyond mass 9, estimates of the triple a reaction 
and a (a n, 7) 9 Be(a, n) 12 C necessary to bridge the mass 
gaps to 12 C and beyond suggest that these three-body 
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Fig. 2. Contours of the final 4 He mass fraction as a func- 
tion of photon to baryon number density ro 7 /nb and co- 
moving dynamical timescale Td yn . 

reactions do not have sufficient time to produce heavy el- 
ements in abundance (see however Section 3). 

The bulk of nucleosynthesis occurs at temperature 
r„ uc ~ 0.1 - 0.2MeV//cB, corresponding to T ~ 10 - 100 
depending on the parameters, and is over by T nuc ~ 
0.06MeV/fce. Nucleosynthesis is favored by a large tdyn 
and a small n 7 /nb, and it is thus more efficient as rj de- 
creases (meaning higher baryon density), as 9 decreases, 
or as r D increases if tdyn oc r Q (since a large expansion 
timescale is more critical to nucleosynthesis than a small 
entropy). Obviously the efficiency of nucleosynthesis is 
maximal for X n /X p — 1 initially, and less efficient as 
this ratio departs from unity. A contour plot of the fi- 
nal mass fraction of 4 He produced is shown in Fig. |^ as a 
function of td yn and n^/n^. In this figure, it was assumed 
that initially X n /X p = 1. It shows that nucleosynthe- 
sis is efficient in most of parameter space: most available 
neutrons and protons are bound into 4 He provided the 
dynamical timescale td yn > 10~ 4 sec if n 7 /nb < 10 4 , or 
idyn > 10~ 3 sec if n-y/nt, < 10 5 . One can understand an- 
alytically these results by comparing the expansion rate 
^dyn t° t ne deuterium formation rate at temperature T 
(nb oc T~ 3 ): i?p( n , 7 )d = (op( n , 7 )dv)A n (l - X n )n h ~ 3 x 
W 4 T-! MeV r;y%^ /2 Lli% ] ) X n (l - X n )sec-\ Indeed 
deuterium production occurs if td yn > 10 4 sec. The de- 
pendence on the entropy comes through the reverse pho- 
todisintegration rate which inhibits the production of deu- 
terium hence of heavier elements. 

One expects nucleosynthesis to be more efficient than 
above if, for instance, magnetic energy represents a signif- 
icant fraction xb of the total energy. The photon number 
density hence n^/n^ are then smaller by (1 — ib) 3 ' 4 while 
the initial temperature is only reduced by (1 — xb) 1 ^ 4 - The 
present calculation also assumes a homogeneous steady 
wind, but inhomogcncitics in the wind, e.g. shells or 
blobs, would increase the efficiency of nucleosynthesis, as 
n 7 /rib oc i? -1 / 4 where E is the local energy content. 



3. Discussion. 

The feedback of nuclear binding energy release on the fire- 
ball evolution has been neglected in the above calcula- 
tions. This is justified since the ratio of entropy density 
released to total entropy density (assuming this latter is 
dominated by the photons) if all neutrons and protons 
are instantaneously converted to 4 He at T rj 0.1 MeV/fcs 
reads dS/S rj 20(fc B T/0.1MeV)- 1 (™7/™b)~\ with dS = 
n b B i /4k B T and B 4 = 28.3 MeV the 4 He binding en- 
ergy. Note that in some rather extreme parts of parameter 
space, the prerequisite conditions for successful r— process 
nucleosynthesis could be satisfied: initial electron fraction 
Y c = X p < 0.5, idyn <; 10~ 3 sec, n 7 /rib ~ 10 2 and initial 
temperature To ~ 1 MeV/fce (e.g. Hoffman et al. 1997; 
Meyer & Brown 1997). This presumably requires a highly 
beamed ejecta 6 < 0.1 from a compact source r < 10 6 cm, 
with a significant fraction of magnetic energy, and/or a 
low limiting Lorentz factor 77. This latter condition can be 
fulfilled in the outer parts of a jet with inhomogeneous 
baryon load as proposed recently (Zhang & Meszaros 
2002). Successful r— process nucleosynthesis might thus be 
able to operate in 7— ray bursts outflows, and further stud- 
ies using dynamical r— process codes appear mandatory. 
Other connections between the site r— process site and 
7— ray bursts have been proposed by Eichler et al. (1989), 
Levinson & Eichler (1993) and Cameron (2001). 

An important consequence of successful nucleosynthe- 
sis to 4 He is to keep neutrons tied to protons and prevent 
their kinematical decoupling when nuclear scattering be- 
comes ineffective. As shown by Bahcall & Meszaros (2000) 
and Meszaros & Rees (2000), such decoupling occurs be- 
fore the Lorentz factor has saturated to its limiting value 77 
provided rj > 480i5o 4 # ~^ 2 r~7^ 4 Y c 1 ^ 4 and leads to relative 
velocities between neutrons and protons ~ c; in turn this 
leads to pion production in inelastic collisions and thus 
to a 5 — 15 GeV neutrino signal. Meszaros & Rees (2000) 
have further shown that transverse diffusion of neutrons 
in inhomogeneously baryon loaded fireballs can lead to an 
appreciable multi-GeV neutrino signal for lower values of 
r\. However these studies assumed that nucleosynthesis did 
not occur. In fact, the inelastic collisions occur in these sce- 
narii at large radii well after nucleosynthesis, hence only 
the surviving free neutrons will be able to decouple. 

Consider first an homogeneous fireball. In the absence 
of nucleosynthesis, the neutrino signal is proportional to 
X p X n = Y e (1 — Yc), and is thus maximal when Y e ~ 0.5. 
The neutron— 4 He collision cross-section is higher than the 
n— p cross-section by the geometrical factor 4 2 / 3 but there 
are only X4/A 4 He nuclei per baryon (X4 denotes as before 
the mass fraction). The neutrino signal produced after nu- 
cleosynthesis (X4 7^ 0) is then a factor r of that produced 
when X4 = (no nucleosynthesis) with: 



2(i-r e ; 



x A 



4Ye (1 - Y c ) 



(2Y C - 0.685A 4 ) 



it was assumed that the outcome of nucleosynthesis is 
only n, p and 4 He. This reduction factor takes the follow- 
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ing values: for Y c ~ 0.5, r ~ 0.33 for A 4 = 0.5, r ~ 0.09 
for X4 = 0.80. If y c > 0.5 and nucleosynthesis is maxi- 
mally efficient, all neutrons are bound into He and obvi- 
ously r = 0. By combining the above relation with Fig. ^| 
it is possible to obtain estimates of r for various 7— ray 
bursts parameters. As an example, r\ = 500 (so that kine- 
matical decoupling of free neutrons occurs before satura- 
tion), L = lQ 50 ergss _1 , 6 = 0.1, r Q = 0.3 x 10 7 cm and 
£dyn = 10 -3 sec gives X4 ~ 0.70 and the neutrino signal is 
reduced by a factor w 6.4. Note that even a small reduc- 
tion factor is significant as the neutrino signal predicted 
is of the order a few events per year for a km 3 neutrino 
telescope and does not exceed a dozen events per year. 

If the jet is inhomogenous, say made of an inner baryon 
poor jet surrounded by a baryon rich outer shell (Eichler 
& Levinson 1999, Meszaros & Rees 2001), the above dis- 
cussion still applies if the outer baryonic wind originates 
from the central engine with a temperature T ~ O(MeV) 
similar to that of the central jet. In effect nucleosynthesis 
is then very efficient in the outer shell since its Lorentz 
factor and entropy are both lower than in the jet. If the 
surrounding shell is "cold" , or if it is neutron rich and the 
inner jet proton rich (or vice- versa) one can circumvent the 
above argument. However in the case of a jet punching its 
way through a collapsar progenitor atmosphere Meszaros 
& Rees (2000) have shown that it leads to a very low neu- 
trino signal. 

Finally, if ultra- high energy cosmic rays are accelerated 
in 7— ray bursts (Levinson & Eichler 1993; Vietri 1995; 
Waxman 1995; Waxman 2001 and references), one expects 
in the present context a significant fraction of these parti- 
cles to be 4 He nuclei. However nuclei are subject to pho- 
todisintegration with cross-section ~a few mb for photon 
energies > 20MeV in the nucleus rest frame, and dur- 
ing acceleration in internal shocks and production of the 
7— ray signal a fraction of them will be disrupted. 

The calculation of the fraction of nuclei dissociated 
as a function of their energy and shock radius is in it- 
self similar to that performed by Waxman & Bahcall 
(1997), Guetta et al. (2001) for the production of a 
100 TeV neutrino signal from pion production of accel- 
erated protons. At the highest energies En c > 10 15 eV 
(observer frame), one can thus write the photodisinte- 
gration rate R = A~ 1 dA/dt ~ U- y cq{a 1 Yi e ) / Ae^, where 
XJ 1 = L 1 /4ir9 2 r 2 ri 2 c is the photon energy density, eb ~ 
1 MeV is the observed break energy in the 7— ray spec- 
trum, and (cr 7 He) is the inverse energy weighted photo- 
disintegration cross-section, (ct 7 h ) — 2.8 mb accounting 
for 1 and 2 nucleon loss with respective branching ratios 
80%, 20%. The optical depth to photodisintegration thus 
reads r w Rt dyn w L 50 e^r^ri^ (e h /1 MeV) -1 , with 
ri4 = r/10 14 cm the radius of emission. 

Since features of temporal width St = 5i _2 x 10~ 2 sec 
are emitted at radius 7-14 ~ 0.65i_2^7foo ; this optical depth 
is unity where features of width St ~ 10 -2 sec are emitted. 
High energy nuclei accelerated in shocks at smaller radii 
are photodisintegrated, while those accelerated at larger 
radii are unharmed. One can also show that low energy 



nuclei are not photodisintegrated even at small radii. The 
overall 4 He nuclei energy spectrum is thus subject to pho- 
todisintegration, acceleration and reacceleration in shocks 
and adiabatic losses. Acceleration at large radii is in any 
case crucial to overcome expansion losses (Waxman 2001). 
It thus seems reasonable to expect that a sizeable frac- 
tion of He nuclei would be present in the escaping ultra- 
high energy radiation and the measurement of the chemi- 
cal composition might provide further information on the 
nuclear processes at work in the 7— ray burst. Detailed 
signatures of these ultra-high energy 4 He nuclei will be 
presented elsewhere. 

Note added: While this paper was being refereed, a sim- 
ilar study by Pruet et al. (2002) appeared, with similar 
conclusions as to the efficiency of nucleosynthesis. 
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